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NNLO BFKL Pomeron eigenvalue in N=4 SYM 


Nikolay Gromov“’^, Fedor Levkovich-Maslyuk“, Grigory Sizov“, 

Mathematics Department, King’s College London, The Strand, London WC2R 2LS, UK 
^ St. Petersburg INP, Gatchina, 188300, St. Petersburg, Russia 

We obtain an analytical expression for the Next-to-Next-to-Leading order of the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) Pomeron eigenvalue in planar A/” = 4 SYM using Quantum Spectral Curve 
(QSC) integrability based method. The result is verified with more than 60 digits precision using 
the numerical method developed by us in a previous paper. As a byproduct we developed a general 
analytic method of solving the QSC perturbatively. 


I. INTRODUCTION 


QCD is notorious for being hard to explore analyti¬ 
cally: perturbative calculations become impossibly com¬ 
plex after first few loop orders. However, there are 
regimes in which one can probe all orders of perturbation 
theory analytically. The Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation is applicable in processes like Deep In¬ 
elastic Scattering or hadronic dijet production, which are 
characterized by a presence of at least two widely sepa¬ 
rated energy scales. The large logarithm of ratio of these 
energy scales Ay enters into perturbative expansion, so 
in order to make sense of the perturbation theory one has 
to resum powers of Ay in every order of perturbation the¬ 
ory. Schematically these large corrections exponentiate 
to an expression of the form 
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where g = with A being’t Hooft coupling constant. In 
particular this can be done in the case of the high-energy 
hadron-hadron scattering [1]. In this case Q absorbs the 
trivial dependence on the transverse momenta and the 
nontrivial part x is the so-called LO BFKL eigenvalue 
[2, 3] 

jn, X ,fn + l + iv\ /n + l — 

X = - - - j-i! - — 

. 

In this paper we focus on the case n = 0. Taking into 
account the Next-to-Leading, Next-to-Next-to-Leading 
corrections we get a similar structure, where the BFKL 
eigenvalue x ii^ th® exponent gets corrected. One usually 
introduces j{iv), related to the BFKL eigenvalue as 
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The Next-to-Leading BFKL was obtained after 9 years 
of laborious calculations in [4-8]; the result in modern 
notation is presented below in the text (3). The correc¬ 
tions turned out to be numerically rather large compared 
to the LO, which makes one question the validity of the 
whole BFKL resummation procedure and its applicabil¬ 
ity for phenomenology. 


This and other indications make it clear that just NLO 
may not be enough to match experimental predictions. 
It is important to understand the general structure of 
BFKL expansion terms and this paper is concerned with 
NNLO BFKL eigenvalue in A/” = 4 SYM — a more sym¬ 
metric analog of QCD. Notably, it was observed in [ 8 ] 
that the Y = 4 SYM reproduces correctly the part of 
the QCD result with maximal transcendentality. In par¬ 
ticular the LO coincides exactly in the two theories. 

Another way of extracting the Pomeron eigenvalue, 
technically more convenient, is due to the observation 
of [ 6 ] who reformulated the problem in terms of a cer¬ 
tain analytical continuation of anomalous dimensions of 
twist-2 operators. Fortunately, in planar Y = 4 SYM 
the problem of computing the anomalous dimensions is 
solved for finite coupling and any operator by the Quan¬ 
tum Spectral Curve (QSC) formalism [9, 10]. 

In order to obtain the BFKL eigenvalue in A/” = 4 
SYM from the anomalous dimension of twist operators 
we consider the dimension A (5) of twist-two operator 
O = TrZD^Z. The inverse function S'(A) is known to 
approach —1 perturbatively for A in the range [—1,1] and 
thus the map to the BFKL regime is given by A = zz/ and 
j = 2-|-S'(A). Then the goal is to compute y(^) as a series 
expansion in g^. Indeed, from the QSC formalism it was 
shown in [ 11 ] that one reproduces correctly the LO (2). 
Here we use some shortcuts to the direct approach of [11] 
to push the calculation to NNLO order, which already 
gives useful new information about the QCD result. 

An essential for us observation was made in [12] [25] 
where it was pointed out that both LO and NLO results 
can be represented as a simple linear combination of the 
nested harmonic sums. Let us stress again that in our 
notation A is the full conformal dimension of the twist- 
two operator, related to the anomalous dimension 7 as 
A = 2 -|- 5 '-|- 7 . Then the expansion of j(A) can be written 
as 
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with the two hrst known orders given by [ 12 ] 
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where 


the form 
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We define harmonic sums for non-integer and negative 
arguments by the standard widely accepted prescription, 
namely analytical continuation from positive even inte¬ 
ger values as in [21, 22]. These analytically continued 
sums, which we denote as S'ai,a2....j are denoted by S~^ in 
[21], see e.g. Eq. (21) in that paper. A compatible but 
more general definition is given in [23]. This prescrip¬ 
tion for analytic continuation is also implemented in the 
Mathematica files attached to the present paper [20]. 

We assume the NNLO order can also be written in this 
form. After that we only have to fix a finite number of 
coefficients which we do by expanding the QSC around 
some values of A where the result simplifies. Then we 
verify our result by comparing it with extremely high 
precision numerical evaluation proving this assumption 
to be correct. 


II. QUANTUM SPECTRAL CURVE 
GENERALITIES 


As it was already mentioned in the introduction there 
is a known relation between the anomalous dimensions of 
the twist-2 operators and the BFKL pomeron eigenvalue. 
Here we describe the Quantum Spectral Curve (QSC) so¬ 
lution of the spectral problem - a simple set of equations 
giving the full spectrum of the anomalous dimensions of 
the theory developed in [9, 10]. Below we limit ourselves 
to the sl{2) sector of the theory. 

The simplest ingredient of the QSC is a set of 4 func¬ 
tions Pa of the spectral parameter u which can be con¬ 
veniently written as a convergent series expansion 

u + y/u— 2gy/u + 2g 
2^ ■ 
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We see that Pa has a branch cut and is power-like at 
infinity. The constants Ma control the global charges of 
the state. For the case of the twist-2 operators, Ma = 
{2,1,0,—!}. The problem of solving the QSC consists 
in finding the coefficients Ca^n- They can be fixed in the 
following steps [13]: 

First, find 4 linear independent analytic in the upper 
half plane solutions of the linear finite difference equation 
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In the generic situation, it is always possible to choose 
the 4 solutions of that equation in this form. In what 
follows we assume this to be done. We require in addition 
that Mi encode the conformal charges of the operators 
i.e. 2M, = (A - S' -f 2, A -k S', -A - S' -k 2,5 - A}, 
where A is the dimension of the operator in question. 
This requirement fixes some of the leading coefficients 
Ca,o- There is an obvious rescaling freedom of Qa\i which 
can be partially fixed by requiring that 

■ (7) 

Next, one finds 4 Q-functions by dividing either side 
of (5) by Pa(u) 

Qj(m) = -P’’(u)Qi,\Au + */2) • (8) 

They are now given implicitly in terms of the coefficients 
Ca,n- The main constraint comes from the condition that 
the analytic continuation which we denote as Qi is a 
linear combination with periodic coefficients of Qi them¬ 
selves 


Qi = uJijQ^ =-'P’’{u)Qb\i{u + i/2) (9) 

where P^(m) is the same as P^(w) with x replaced by 1/x. 
In particular from (9) we have 

Qi(w) = Wi2Q3(m) + Wi4Qi(m) — a;i3Q2(u) , (10) 

Q3(m) = W34Qi(m) - Wi4Q3(m) + Wi3Q4(u) . 

As we will see, uJij can be eliminated. To show this we will 
need only to know that ojij is Aperiodic, anti-symmetric 
and should satisfy 

LOijUjy^ = S/ , W23 = Wi4 . (II) 

For all physical operators ujij should go to a constant at 
large u, however it was emphasized in particular in [13] 
and [14] that for non-integer S one should allow for an 
exponential growth of a;24 as otherwise the system has no 
solution. Note that (11) implies that = ojis should 
decay exponentially at infinity. It is also known that 
a;i4 decays at infinity [13]. Condition (9) in fact imposes 
infinitely many constraints on the coefficients Ca,n fixing 
them completely as well as the function A(S) or iS'(A). 


III. ANALYTICAL DATA FROM QSC 


Qa|i(w+§)-2a|i(u-|) = -Pa(u)PAu)Qbli(u+^) (5) 

where i labels the 4 solutions. Here and everywhere in 
this paper indices are raised with a 4 x 4 matrix = 
(—I)“da.5-&- The solutions can be always chosen to have 
“pure” asymptotics, which means that with exponential 
precision the large u (asymptotic) expansion of Qa\i has 


We describe now the details of our analytical method. 
We will focus on some particular points Aq = 1,3, 5,7. 
It can be seen already from the LO (2) that the function 
^(A) is singular at these points, however the coefficients 
of the expansion are relatively simple and are given by 
^-functions. We will perform a double expansion first in 
g up to the order and then in d = A — Aq. 
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a. General iterative procedure for solving QSC. We 
describe a procedure which for some given (or equiva¬ 
lently Ca,n) takes as an input some approximate solution 
of (5) valid up to the order e" (where e is some 
small expansion parameter) and produces as an output 
new Q^|- accurate to the order The method is very 
general and in particular is suitable for perturbative ex¬ 
pansion around any background. 

Let dS be the mismatch in the equation (5), i.e. 

(12) 

where dSa\i is small ~ e". We can always represent the 
exact solution in the form 

Qa\^{u) = Qa\l(^) + Wiu + §) 2l|’](w) (13) 

where the unknown functions &/ are also small. After 
plugging this ansatz into the equation (12) we get 

(b-\u) - bf{u + i)) = dSa\^ + dSaijbf . (14) 

Since bf is small it can be neglected in the r.h.s. where 
it multiplies another small quantity. Finally multiplying 
the equation by and using (7) we arrive at 

bliu + *) - (u) = (u + §) + 0(e"") • 

We see that the r.h.s. contains only the known functions 
dS and and does not contain b which means that the 
original 4th order finite difference equation is reduced 
to a set of independent 1st order equations! In most 
interesting cases the first order equation can be easily 
solved. After is found one can use (8) to find Q^. 

b. Iterations at weak coupling. For our particular 
problem we will take either e = g or e = S. Applying 
this procedure a few times we generate Qi for sufficiently 
high order both in g and in 5. Finally, by “gluing” 
and Qi on the cut we hnd Ca,n and S{A) also as a double 
expansion. 

For the above procedure we need the leading order 
expect that to the leading order in g the 
solution should be very simple - indeed the branch cuts 
collapse to a point making most of the functions polyno¬ 
mial or having very simple singular structure. Also one 
can use that to the leading order in g functions Pq are 
very simple and are already known from [11] for any A. 
By making a simple ansatz for Qi we found for Aq = 1 
to the leading order 

Qi ~ u, Q 2 ~ 1/m, Q 3 ^1, Q4 ^ 1/m^ . (15) 

For Aq = 3,5,... the solution involves also the r/- 
functions introduced in the QSC context in [15, 16] 

r]si,...,sk{u) ^ _ _ (-y _|_ ■ 

ni>n2...nfc>0 

( 16 ) 


which are related in a simple way to the nested harmonic 
sums. For A = 3 we found 

Ql ~ M^, Q 2 ~ M^7?i, 3 - z - (17) 

Zu 

2 1 2 i 1 

Qs-U 771,2 -lU- Q 4 ^ M 771,4 - - - ^ , 

which reflects the general structure of the expansion of 
Qi around integer A’s which contain only 771 , 2 , 771,3 and 
771,4 with polynomial coefficients. As it was explained in 
[15, 16] the 77 -functions are closed under all essential for 
us operations: the product of any two 77 -functions can be 
written as a sum of 77 -functions, and most importantly 
one can easily solve equations of the type 

/(M-f z) -/(m) = m”77^,....,s, (18) 

for any integer n again in terms of a sum of powers of 
u multiplying 77 -functions (which we call 77 -polynomials). 
For example for n = — 1 and /c = 1, si = 1 we get 
/ = —772 — 771,1 etc. Thus for these starting points we 
are guaranteed to get 77 -polynomials on each step of the 
general procedure described above. 

Proceeding in this way we computed Qi up to the order 
g^ and for A = 3, 5, 7. After that we fix the coeffi¬ 
cients in the ansatz for from analyticity requirements 
described below. 

c. Fixing remaining freedom. Here we will describe 
how to use Qi found before to finally extract relation be¬ 
tween S and A and the constants Ca,n- This is done by 
using a relation between Qi and their analytical contin¬ 
uations Qi. On the one hand we have the relation (10). 
On the other hand we can use the u —^ —u symmetry[26] 
of the twist-2 operators to notice that Qi(— u) should 
satisfy the same finite difference equation as Qi(M) and 
thus we should have Qi(M) = 0 /(M)Qj(— m) where flj(u) 
is a set of periodic coefficients. As Qi(M) has a power-like 
behavior at infinity, flj (u) should not grow faster than a 
constant. Furthermore, since Qi has a definite asymp¬ 
totic ( 6 ) only diagonal elements of 0*(m) can be nonzero 
at infinity. Combining these relations we find 

Q 2 i(m) = a)iQi(-u) , A = 1,3, (19) 

where a\ = uJAiX'‘^^k Aperiodic (as a combination 
of Aperiodic functions), analytic (as both Qa{u) and 
Qa(— m) should be analytic in the lower-half-plane) and 
growing not faster than a constant at infinity which im¬ 
plies that they are constants. Furthermore most of them 
are zero because only a;i 2 , 0734 and fl® are non-zero at 
infinity. Thus we simply get 

Qi(m) = q;i 3 Q 3 (-m) , Q3(m) = Q!3iQi(-m) • (20) 

Next we note that if we analytically continue this relation 
and change m —>■ — u we should get an inverse transfor¬ 
mation which implies 013 = 1 / 0:31 = a. The coefficient 
a depends on relative normalization of Qi and Q 3 . Let 
us see how to use the identity ( 20 ) to constrain the con¬ 
stants Ca,i. We observed that all the constants are fixed 
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from the requirement of regularity at the origin of the 
combinations Qi + Qi and , which now can be 

■y/u2-4g2 


written as 


Qi(m) + aQ 3 (-u) = reg 


Qi(u) - a;Q3(-u) 
\/u^ — 


= reg . 


J 


This relation is used in the following way: one first ex¬ 
pands in g the l.h.s. and then in u around the origin. 
Then requiring the absence of the negative powers will 
fix a, all the coefficients Ca,n, and the function A(5)! So 
we can completely ignore Q 2 , and Q 4 in this cal¬ 
culation. This observation can be used in more general 
situations and allows avoiding construction of and 
in particular can simplify the numerical algorithm of [13] 
considerably. 


d. Constraints from poles. We use the procedure described above to compute the expansion of S'(A) around 
Aq = 3, 5, 7. In particular for A = 5 -I- e we computed the first 8 terms 


X 


NNLO 


1024 64 (47r2 - 33) 16 (-36C3 + + 3l) -288C3 + - 1671^ - 296 


3e3 


( 21 ) 


[20 (47r2 - 75) Cs + 6300C5 + - 2157r2 -p 285] -b ... . 


The terms with e, e^, and which we also evaluated explicitly are omitted for the sake of brevity. We also reproduced 
expansions extracted from [17] for A = 1. In our calculations we used several Mathematica packages for manipulating 
harmonic sums and multiple zeta values [18]. 


IV. THE RESULT 


By observing (4) for LO and NLO we notice that 
the transcendentality of these expressions is uniform 
if one assigns to Sai,...,ak transcendentality equal to 

fc 

^ ja^-j. The principal assumption of our calcula- 
1=1 

tion states that ^ 3 ( 0 ;) can also be written as a lin¬ 
ear combination of nested harmonic sums with coef¬ 
ficients made out of several transcendental constants 
TT^,log( 2 ),^ 3 ,^ 5 ,Li 4 (i) ,Li 5 (l) of uniform transcenden¬ 
tality 5. The final basis obtained after taking into ac¬ 
count the constants contains 288 elements. 


have a unique solution presented below: 


256 8 2 2 2 

S-4.S1 S-3S2 383^-2 35'_3,i4 

^ + — -^ ^- 2 . 1,1 

35 _ 2^3 ^5 , 5 - 2^1 ^2 

+ <52,-2,! + 1><5-2, 1 , 1,1 -- 5 - H-- 

+C 3 


( 22 ) 


'5-2,1 75-3 

5-2 5i 5 i52' 

TT^ 

'25-1 

5i1 

00 

00 

12 48 . 


45 

96 


75-1,1 78.2 7S.iSi 82 

4 8 4 16 


4 

.4 on 


2Li4(i) - 

12 ^ 12 


(5-1 - 5i) 


log_^ _ 7 r 2 log32 _ 2Tr'^ log 2 _ 49(5 _ m 

36 45 24 32 ^421 


60 


The simplicity of the final result is quite astonishing: only 
37 coefficients out of 288 turned out to be nonzero. Fur¬ 
thermore, they are signihcantly simpler than the coeffi¬ 
cients appearing in the series expansion around the poles 
(21). These are all clear and expected indications of the 
correct result similar to what was observed in the usual 
perturbation theory [19]. In addition we also performed 
the numerical test described below. 


Hence we build the linear combination of these basis 
elements with free coefficients and constrained them by 
imposing the expansion at A = 1,3, 5,7 to match the 
results of the analytic expansion of QSC (in particular, 
requiring (21)). This gave an overdefined system of linear 
equations for the unknown coefficients which happen to 


V. NUMERICAL TESTS 

Using the method of [13] we evaluated 40 values of 
spin 8 for various values of the coupling g in the range 
(0.01,0.025) with exceptionally high 80 digits precision 
and then fit this data to get the following prediction for 
the N"LO BFKL coefficients at the fixed value of A = 
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0.45: 



value 

error 

N^LO 

10774.6358188471766379575931271924 

56995929170948057653783424533229 

10-61 

N^LO 

-366393.20520539170389379035074785 

44549935531959333919163403836 

10-56 

N^LO 

1.33273635568112691569404431036982 
8561521940588979476878854 x 10^ 

io-®i 

N®LO 

-4.9217401366579165009139555520750 
70060721450958436559876 x 10® 

10-47 


We found that our result (22) reproduces perfectly the 
first line in the table within the numerical error 10“®^ 
which leaves no room for doubt in the validity of our 
result. 


VI. SUMMARY 

In this letter we have applied the Quantum Spectral 
Curve method [9, 10] to the calculation of the NNLO 
correction to the BFKL eigenvalue. We check our re¬ 
sult numerically with high precision using the algorithm 
developed in [13] and gave numerical predictions for a 
few next orders. We also developed a general efficient 
analytic method suitable for systematic perturbative so¬ 
lution of QSC. 


There are numerous packages such as [18] available for 
the evaluation of the nested harmonic sums. Yet to sim¬ 
plify future applications of our results we attached a small 
Mathematica notebook [20] which allows to evaluate our 
result numerically and also generate expansions about 
singularities and at infinity. 

We hope that our findings could shed some light on the 
QCD counterpart of our result and resolve some myster¬ 
ies shrouding the BFKL physics. 
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